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Preface 

This  thesis  represents  six  months  of  research  in  the  Particles 
and  Fields  Division  of  the  Space  Physics  Laboratory  at  Aerospace 
Corporation  as  a  part  of  the  AFIT  masters  degree  program.  During 
this  time  I  analysed  satellite  data  to  gain  information  about  outer 
belt  electron  behavior  and  the  magnetosphere. 

Knowledge  of  the  radiation  belt  is  important  both  to  science  and 
die  Air  Force.  To  science  it  represents  a  better  understanding  of  our 
environment.  To  the  Air  Force*  which  is  presently  concerned  with 
manned  orbital  flight  it  represents  knowledge  about  one  of  die  hazards 
encountered  in  space.  It  is  important  to  precisely  know  the  radiation 
hasards  in  order  to  protect  both  men  and  equipment  from  damage. 

The  personnel  of  the  Space  Physics  X*ab  provided  me  with  im¬ 
mense  support  and  for  this  I  am  deeply  indebted.  Virtually  every  in¬ 
dividual  played  some  part  in  making  this  thesis  possible  and  I  wish  to 
thank  them  all.  I  would  specifically  like  to  recognize  Dr.  Albert 
Vampola  and  Dr.  George  Paulikas.  They  not  only  gave  me  key  advice 
and  guidance  throughout  my  research  but  made  my  stay  a  learning  ex¬ 
perience  extending  well  beyond  the  area  of  my  specific  topic.  I  thank 
them  both  for  their  valuable  time  and  effort. 

I  would  like  to  specially  thank  Lt.  Col.  Donald  L.  Evans  both  for 
his  valuable  assistance  as  my  AFIT  thesis  advisor  and  for  making 
this  valuable  learning  experience  possible.  I  also  acknowledge  the  as* 
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sistance  extended  me  by  Lt.  David  J.  Evans  of  SA&SO* 
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Abstract 

The  high  latitude  electron  cutoffs  are  determined  from  electron 
flux  data  collected  by  the  magnetic  spectrometer  aboard  the  polar  or¬ 
biting  OV3-3  satellite  (1966-70A).  The  cutoffs  are  examined  as  a 
function  of  nine  energies  from  300  k.eV  to  2.3  MeV  and  are  found  to  be 
a  function  of  energy  for  all  local  times  and  to  exhibit  an  energy  depen¬ 
dent  diurnal  variation.  This  phenomenon  is  explained  with  the  use  of 
the  models  of  Taylor  and  Hones  and  of  Williams  and  Mead.  The  lower 
energy  cutoffs  occurring  at  higher  latitudes  is  purported  to  result 
from  kinks  in  the  nightside  magnetic  field  lines  and  the  limiting  ef¬ 
fect  of  rigidity.  The  greater  diurnal  variation  for  the  lower  energy 
cutoffs  is  thought  to  be  the  result  of  field  distortion  with  the  different 
energies  labeling  different  L  shells.  The  cutoffs  show  no  altitude  de¬ 
pendence  and  no  well  defined  relationship  with  magnetic  activity  as 
represented  by  either  Dst  or  Kp. 
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ENERGY  DEPENDENCE  OF 
HIGH  LATITUDE  ELECTRON  CUTOFFS 

I.  Introduction 

The  radiation  belt  surrounding  die  earth  was  first  discovered  by 
Van  Allen,  Ludwig,  Ray  and  Mcllwain  in  1953  (Ref.  3:220).  The  belt 
was  shown  to  contain  both  electrons  and  protons  (Ref.  4:16).  A  map¬ 
ping  of  die  electrons  in  die  radiation  belt  show  that  for  high  energy 
electrons  the  belt  consists  of  two  zones,  while  for  lower  energies  die 
distribution  is  almost  continuous  to  about  10  earth  radii  at  die  equator. 
For  energetic  electrons  the  inner  zone  extends  to  around  2  earth  radii 
and  the  outer  zone  from  about  3  to  10  earth  radii.  This  is  shown  in 
Figure  1.  In  the  inner  and  outer  zones,  source,  loss,  acceleration 
and  diffusion  mechanisms  still  remain  quantitatively  uncertain. 

The  present  study  is  an  investigation  of  trapped  electron  beha¬ 
vior  in  the  outer  fringes  of  the  outer  zone  at  low  altitudes.  The  high 
latitude  cutoff  of  trapped  electrons  is  examined  to  gain  information 
about  outer  zone  electron  behavior  and  the  magnetosphere.  Observing 
the  outer  zone  electrons  yields  clues  as  to  what  is  going  on  in  the 
outer  portions  of  the  magnetosphere.  The  high  latitude  cutoffs  are  im¬ 
portant  in  examining  the  trapping  boundary  or  the  limit  of  closed  field 
lines  in  the  magnetosphere.  Figure  2  shows  the  relationship  between 
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the  high  latitude  cutoff  and  the  field  lines  of  the  magnetosphere.  At 
present*  the  outer  part  of  the  magnetosphere  is  die  least  understood. 
Most  magnetospheric  models  are  valid  close  in  to  the  earth,  but  this 
validity  decreases  rapidly  beyond  several  earth  radii.  Information  a- 
bout  die  limit  of  closed  field  lines  can  be  used  as  a  boundary  condition 
to  evaluate  and  improve  magnetospheric  models* 

This  investigation  is  a  systematic  examination  ol'  the  high  lati¬ 
tude  cutoffs  from  approximately  1 50  orbits  of  the  OV3-3  polar  orbiting 
satellite  (1966-70A).  The  cutoffs  are  analyzed  as  a  function  of  nine 
energies  from  300  keV  to  2.3  MeV  and  relationships  with  magnetic  ac¬ 
tivity  are  examined.  The  computer  is  used  extensively  to  facilitate 
examination  of  the  available  data. 
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II.  Background 

Particle  motion  in  the  magnetosphere  is  briefly  reviewed  as  a 
means  of  clarifying  the  significance  of  the  high  latitude  cutoff  of  elec¬ 
trons.  The  magnetic  field  of  the  earth  can  be  approximately  repre¬ 
sented  as  a  dipole,  as  shown  in  Figure  3.  A  charged  particle  moving 
in  such  a  field  will  have  three  distinct  motions;  the  particle  will  spi¬ 
ral  about  a  field  line  (cyclotron  motion),  bounce  back  and  forth  along 
the  field  line  between  two  mirror  points,  and  drift  around  the  earth. 
(Ref.  4:23).  These  three  motions  are  depicted  in  Figure  4.  The  ap¬ 
proximate  periods  associated  with  these  motions  for  a  1  MeV  electron 
with  a  2000  km  altitude  at  the  equator  are  10"^  sec  for  the  cyclotron 
motion,  10" 1  sec  for  the  bounce  motion  and  60  min  for  the  drift  mo¬ 
tion  (Ref.  4:34).  This  close  association  between  these  trapped  par¬ 
ticles  and  the  magnetosphere  make  the  behavior  of  the  particles  use¬ 
ful  in  investigating  the  magnetosphere. 

The  assumptions  that  the  magnetic  field  is  slowly  varying  with 
time  and  space  yields  three  adiabatic  invariants.  The  magnetic  mo¬ 
ment,  the  length  of  a  particle's  bounce  path  along  a  field  line,  and  the 
magnetic  flux  enclosed  by  the  drift  path  can  be  assumed  constant  if  the 
time  scale  of  any  variation  and  spatial  scale  of  any  inhomogeneity  in 
the  magnetic  field  are  large  compared  to  those  characteristic  of  the 
cyclotron,  bounce,  and  drift  motions,  respectively.  The  result  of  the 
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magnetic  moment  being  constant  is  to  require  that  the  flux  through  a 
cyclotron  orbit  remain  constant  and  that  a  specific  particle  will  al  - 
ways  mirror  at  the  same  value  of  the  magnetic  field. 

These  three  anabatic  invariants  are  three  constants  of  motion 
that  one  can  use  to  describe  a  particle's  motion  in  the  actual  dynamic 
field  of  the  earth,  assuming  B  changes  little  with  time  and  space. 
However,  B  does  change  with  time  and  space,  and  on  occasion  one  or 
more  of  the  invariants  may  be  violated.  The  third  is  the  most  likely 
of  the  three  to  be  violated  since  changes  of  the  order  of  an  hour  are 
not  unusual  during  magnetically  active  periods.  There  is  evidence 
that  the  other  invariants  may  also  be  violated. 

To  simplify  the  discussion  and  representation  of  trapped  par¬ 
ticles  in  the  imperfect  dipole  field  of  the  earth,  Mcflwain  proposed 
using  the  value  of  the  magnetic  field,  B,  and  the  magnetic  shell  para¬ 
meter,  L,  as  coordinates.  The  advantage  of  this  system  is  that  a  set 
of  values  of  B  and  L  correspond  uniquely  to  one  set  of  values  of  B 
and  I,  the  integral  invariant.  A  trapped  particle  ideally  will  always 
have  the  same  L  coordinate  in  this  system  as  it  drifts  around  the  earth. 
For  a  dipole  field,  L  is  equal  to  the  equatorial  distance  from  the  cen¬ 
ter  of  the  earth  to  the  field  line,  and  for  a  nondipole  field,  L  is  the 
generalization  of  the  equatorial  distance,  which  turns  out  to  be  very 
nearly  constant  along  a  line  of  force.  I*  varies  by  leas  than  1  percent 
along  most  field  lines  (Ref.  4:57-58).  Figure  5  shows  how  the  mag¬ 
netosphere  is  described  in  the  B-L  coordinate  system. 
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Associated  with  the  B-L  coordinate  system  is  the  invariant  lat¬ 
itude. 


A  = 


-i  ,1.1/2 

cos  (-) 


0) 


This  latitude  can  be  considered  very  approximately  as  the  latitude 
where  the  L  shell  intersects  the  earth.  For  an  infinite  L  value,  A  = 

90  deg;  for  L  =  1,  A  =  0  (Ref.  4:59). 

The  interaction  of  the  solar  wind  and  the  geomagne'dc  field  has 
been  described  by  William  and  Mead.  The  solar  wind  compresses  the 
field  lines  on  the  day  side  of  the  earth  and  extends  those  on  the  night 
side,  as  shown  in  Figure  6.  The  boundary  between  the  solar  wind  and 
the  magnetosphere  on  the  day  side  is  located  at  about  10  Re  (Ref.  16: 
3024).  On  the  night  side,  the  effect  of  the  solar  wind  is  to  cause  field 
lines  beyond  about  7  Re  to  be  open  instead  of  closed.  Also,  for  a  given 
altitude  above  the  earth,  the  value  of  B  is  greater  on  the  day  side  than 
on  the  night  side  as  a  result  of  the  solar  wind  compression  of  the 
field  lines  (Ref.  16:3025). 

The  high  latitude  cutoff  of  electrons  is  the  invariant  latitude 
above  which  trapped  electrons  are  not  observed.  If  closed  field  lines 
extend  only  to  7  Rfi  on  the  night  side  of  the  earth  and  the  boundary  of 
the  magnetosphere  is  at  10  Re  on  the  day  side,  no  particles  should  be 
trapped  for  L  >  7  on  the  night  side  or  L  >  10  on  the  day  side.  From 
eq.  (1)  we  can  calculate  the  invariant  latitudes  corresponding  to  these 
two  L  values  and  the  theoretical  day  side  and  night  side  high  latitude 
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cutoffs.  They  are  approximately  71.5  deg  invariant  latitude  on  the 
day  side  and  68  deg  on  the  night  side.  Actual  measurements  of  the 
cutoffs,  however,  have  shown  that  these  calculations  are  not  born  out 
by  experiment. 

O'Brien  discovered  in  1963  that  the  high  latitude  cutoff  for  elec¬ 
trons  with  energy  E>40  keV  varied  by  several  degrees  according  to 
the  local  time,  with  the  highest  values  occurring  in  the  local  day.  The 
cutoffs  occurred  at  »  75  deg  around  local  noon  and  a*  69  deg  around 
local  midnight  (Ref.  9:989). 

In  1964,  McDiarmid  and  Burrows,  in  examining  Aloutte  data  for 
electron  energy  E  >40  keV,  found  substantially  the  same  diurnal  vari¬ 
ation  in  the  cutoffs  as  O’Brien.  They  noted  that  the  cutoffs  were  ap¬ 
proximately  symmetrical  about  the  noon-midnight  meridian  (Ref.  6: 
616). 

In  1965,  Williams  and  Mead  examined  the  high  latitude  boun¬ 
daries  for  electrons  of  energies  E  >280  keV  and  E  >1.2  MeV.  They 
found  a  noon-midnight  asymmetry  of  2  l/2  deg  which  corresponds 
closely  with  the  asymmetry  predicted  by  invariant  particle  motion  in 
a  distorted  magnetosphere.  They  concluded  that  the  boundary  was  the 
same  for  electrons  with  energies  E  >40  keV,  E  >280  keV,  and  E  > 

1.2  MeV  at  local  midnight.  However,  the  dayside  cutoff  for  electrons 
E  >40  keV  consistently  extended  to  higher  latitudes  than  those  fox 
electrons  with  energies  E  >280  keV  and  E>1.2  MeV  (Ref.  16:3028). 

Also  in  1965,  Armstrong  examined  cutoffs  from  Injun  3  satel- 
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lite  data  and  found  the  diurnal  variation  to  be  a  function  of  energy.  He 
investigated  electrons  with  energies  E  >40  keV*  E  >230  keV,  and  E  > 
1.6  MeV  and  noted  that  the  diurnal  variation  decreased  with  increasing 
energy  (Ref.  1:2101). 

In  I960,  McDiarmid  and  Burrows  again  examined  the  high  lati¬ 
tude  cutoffs  as  a  function  of  energy  and  found  that  the  E>3.9  MeV  cut¬ 
offs  were  approximately  symmetrical  with  respect  to  the  1100-2300 
hr  meridian  and  exhibited  a  diurnal  variation  of  3  deg.  For  E>35  keV 
electrons  a  diurnal  variation  of  6  deg  was  found,  and  the  cutoffs  were 
not  symmetrical  about  any  meridian.  However,  an  appreciable  dawn- 
dusk  asymmetry  was  evident  (Ref.  7:49).  These  trends  can  be  seen  in 
Figure  7. 

The  role  of  geomagnetic  activity  in  influencing  the  high  latitude 
cutoffs  of  electrons  was  examined  by  Maehlum  and  O’Brien  in  1963. 
They  found  that  during  times  of  enhanced  magnetic  activity  the  boun¬ 
dary  for  electrons  E>40  keV  moved  to  lower  latitudes  for  both  night 
and  day  cutoffs.  The  cutoffs  also  became  quite  sharp  during  storms 
(Ref.  5:977-1000). 

The  present  study  is  concerned  with  clarifying  the  above  obser¬ 
vations.  The  data  used  consist  of  magnetic  spectrometer  measure¬ 
ments  of  Jl  ,  the  electron  flux  perpendicular  to  the  field  line  (<*=  90± 

8  deg),  for  nine  energies  from  300  keV  to  2.3  MeV.  The  instrument 
is  shown  in  Figure  8.  Its  design,  construction,  and  calibration  were 
performed  by  A.  L.  Vampola  at  Aerospace  Corporation  and  were  not 
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a  part  of  this  study.  Details  of  the  instrument's  operation  are  dis¬ 
cussed  in  Appendix  A.  The  vehicle  carrying  the  instrument  was  the 
OV3-3  pclsr  satellite  launched  on  4  August  1966.  The  satellite  had  a 
period  of  136.6  min*  an  inclination  of  81.6  deg*  a  perigee  of  200  nm, 
and  an  apogee  of  2780  nm.  A  schematic  of  the  satellite's  trajectory  is 
shown  in  Figure  9.  The  data  was  merged  with  the  ephemeris  by 
computer*  and  plots  of  versus  universal  time  obtained  for  about  1 50 
orbits.  An  example  of  these  plots  appears  in  Figure  10.  The  high  la¬ 
titude  cutoffs  used  in  this  study  were  obtained  from  these  plots. 
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FIG.  9.  FLIGHT  CHARACTERISTICS  OF  THE  OV3-3  SATELLITE 
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III.  Data  Analysis 

The  problem  of  defining  exactly  where  the  high  latitude  cutoff 
occurs  is  not  new.  McDiarmid  and  Burrows  dealt  with  this  problem 
by  defining  and  using  several  cutoffs.  For  example*  they  defined  a 
background  cutoff  where  the  flux  falls  to  the  background  level*  a 
smooth  cutoff  where  the  flux  profile  shows  a  marked  change  in  charac¬ 
ter*  and  a  sharp  cutoff  where  the  flux  decreases  by  a  factor  of  ten  or 
more  in  a  latitude  interval  of  less  than  half  a  degree  (Ref.  7:50).  Of 
the  three*  the  only  definition  practical  for  this  study  was  the  first,  or 
background,  cutoff.  Defining  the  cutoff  at  the  point  where  the  flux  falls 
to  a  certain  percentage  of  the  maximum  was  considered,  but  any  rea¬ 
sonable  percentage  eliminated  too  much  data.  The  maximum  flux 
values  for  the  highest  energy  channels  were  often  so  low  that  taking 
five  percent  of  these  values  yielded  a  flux  value  below  background* 
making  a  cutoff  point  impossible  to  define.  Taking  ten  or  more  per¬ 
cent  of  the  maximum  made  the  cutoff  point  too  dependent  upon  the  flux 
profile  and  therefore  meaningless.  Thus  the  cutoffs  used  in  this  study 
were  defined  as  the  invariant  latitude  at  which  ,  the  flux  perpen¬ 
dicular  to  the  local  field  line  (ot=  90±_8  deg),  assumed  its  background 
level. 

Using  a  background  cutoff  was  not  without  problems.  The  first 
was  with  cutoffs  that  were  not  extremely  sharp.  Where  the  flux  curve 
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slowly  approached  background,  the  statistical  variation  associated 
with  tiie  low  background  counting  rates  made  determination  of  the 
exact  cutoff  point  difficult.  This  problem  was  overcome  by  extra¬ 
polation  of  the  flux  curve  down  from  an  order  of  magnitude  above 
the  background  level. 

The  second  problem  was  the  effect  of  solar  flare  electrons  (Ref 
14:1254)  on  the  polar  background  levels.  They  increased  the  normal 
background  levels  by  orders  of  magnitude  for  the  lower  energy  chan¬ 
nels.  while  the  higher  energies  were  relatively  unaffected.  This  was 
overcome  by  simple  extrapolation  down  to  the  observed  background 
levels  for  each  channel  and  accepting  the  decreased  reliability  of  the 
data  from  the  affected  orbits.  As  shown  in  Figure  11,  around  a  univer¬ 
sal  time  of  56450  sec  the  shape  of  the  flux  curve  has  been  altered,  and 
it  is  difficult  to  say  exactly  where  trapping  stops. 

Approximately  2800  cutoffs  were  defined.  Each  cutoff,  along 
with  the  associated  parameters  of  interest  from  the  ephemeris,  was 
put  on  cards  to  facilitate  computer  processing. 

The  cutoffs  were  usually  sharper  and  better  defined  during  local 
night  than  during  local  day.  Also,  in  several  instances  low  energy 
(  »  300  keV)  spikes  were  seen  above  the  cutoff  latitude  during  mag¬ 
netically  disturbed  periods.  Even  in  some  cases  where  no  distinct 
spike  was  present,  the  cutoff  shape  was  altered  in  such  a  way  as  to 
suggest  a  spike  being  barely  enclosed  by  the  cutoff.  For  example,  the 
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cutoff  might  appear  pushed  out  and  deformed  just  before  the  flux  fell 
to  background.  Figure  12  shows  these  phenomena.  The  local  times 
during  which  spikes  were  observed  are  listed  in  Table  I.  Most  de¬ 
finite  spikes  occur  around  local  midnight,  while  evidence  of  less  de¬ 
finite  spikes  was  observed  mostly  in  the  morning  before  noon.  No 
conclusion  with  regard  to  the  spikes  was  reached,  and  too  few  were 
observed  to  facilitate  a  detailed  investigation. 

There  is  one  problem  inherent  in  looking  at  cutoffs.  When  the 
cutoff  is  observed  to  be  different  for  any  two  revolutions,  the  obser¬ 
ver  cannot  say  whether  pariicles  moved  to  different  L  shells  or  the 
L  shells  moved  to  different  locations,  because  to  get  values  for  L 
and  A  we  use  a  computerized  magnetospheric  model,  which  does  not 
take  into  account  all  such  variations  in  the  magnetosphere. 

The  character  of  the  data  virtually  dictated  a  statistical  approach. 
The  data  was  taken  at  different  altitudes,  local  times,  and  conditions 
of  magnetic  activity.  In  addition,  coverage  was  such  that  in  few  cases 
are  there  two  consecutive  revolutions  of  data,  and  even  two  cutoffs 
taken  under  similar  conditions  are  few.  In  most  cases  the  data  was 
separated  by  days  and  conditions  varied  accordingly.  It  was  thus  im¬ 
possible  to  relate  the  changes  of  cutoffs  over  an  extended  period  of 
time  to  the  changing  conditions. 

The  local  time  referred  to  in  this  study  is  the  geographic  local 
time  of  the  point  where  the  cutoff  was  observed.  The  expression  used 
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TABLE  I:  Data  on  High  Latitude  Spikes 


Day  of 
Observation 

Local 

Time  (hours) 

Type  of 
Observation 

235 

1415 

definite  spike 

249 

2128 

ft  ft 

250 

0808 

evidence  of  a  spike 

251 

1049 

ft  ft  ft 

251 

1013 

ft  ff  ft 

251 

2246 

definite  spike 

252 

1304 

evidence  of  a  spike 

252 

2105 

definite  spike 

252 

0731 

evidence  of  a  spike 

253 

0730 

ft  ft  ft 

259 

0944 

ft  ft  ft 

271 

2214 

ft  ft  ft 

271 

0551 

ft  ft  ft 

271 

0014 

definite  spike 

273 

0605 

evidence  of  a  spike 

273 

0645 

ft  ft  ft 

277 

2250 

ft  ft  ft 

297 

0930 

definite  spike 

339 

0145 

ft  ft 

339 

0546 

evidence  of  a  spike 
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in  its  calculation  is 

LT  =  UT/3600  +  E/15  (2) 

where  LT  is  the  local  time  in  hours,  UT  is  the  universal  time  in  sec¬ 
onds  and  E  is  the  east  longitude  in  degrees.  Probably  a  more  correct 
parameter  would  have  been  magnetic  local  time.  However,  due  to  the 
time  and  effort  required  to  obtain  this  parameter,  geographic  local 
time  was  used  instead.  McDiarmid  (private  communication)  stated 
that,  in  his  study  of  cutoffs,  he  found  little  difference  between  magne¬ 
tic  and  geographic  local  times. 

From  the  data  on  cards,  initial  plots  were  made  of  cutoffs  ver¬ 
sus  altitude,  local  time,  and  the  3  hr  averages.  Each  energy  for 
each  pole  was  plotted  separately.  Although  there  was  a  good  deal  of 
scattering  of  the  points,  it  was  clear  that  for  the  statistical  approach 
being  used  the  North  Pole  and  South  Pole  data  could  be  superimposed. 
This  was  not  unexpected  since  the  magnetospheric  model  used  by  the 
computer  for  data  reduction  corrected  for  the  north-south  asymmetry 
of  the  magnetosphere. 

This  superposition  yielded  about  double  the  density  of  data 
points  to  work  with.  Running  these  plots  off  again  for  each  energy 
still  left  the  scattering  problem  apparent  in  the  initial  plots.  Exam¬ 
ples  of  the  three  types  of  plots  appear  as  Figures  13  through  15.  Since 
variations  in  magnetic  activity  should  be  approximately  random,  rela¬ 
tions  to  remove  the  altitude  and  local  time  dependences  were  sought. 
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FIG.  13.  PLOT  OF  CUTOFF  versus  ALTITUDE  FOR  CHANNEL  7 
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An  examination  of  the  nine  plots  of  A  versus  altitude  indicated  that 
the  cutoffs  were  independent  of  altitude  in  this  experiment.  In  none  of 
the  nine  cases  did  any  trend  seem  to  be  indicated.  If  any  trend  were 
really  present,  it  should  at  least  be  detectable  in  these  plots. 

The  initial  plots  made  it  clear  that  there  was  a  local  time  de¬ 
pendence.  Half-hour  averages  and  the  associated  standard  deviations 
were  computed  and  plotted.  This  was  done  to  eliminate  any  confusion 
caused  by  the  varying  density  of  points  over  local  time  and  to  prevent 
biasing  of  any  attempted  fits  to  the  data.  As  shown  in  Figure  16,  the 
scattering  still  existed,  although  a  definite  trend  was  visible.  A 
fourth  order  fit  to  this  average  data  failed  to  merge  smoothly  around 
the  endpoints.  Thus,  the  average  values  from  0  to  24  hr  were  repro¬ 
duced  from  24  to  48  hr  and  a  ninth  order  fit  was  found  satisfactory, 
as  shown  in  Figure  17.  Then  a  fourth  order  fit  to  this  ninth  order 
curve  over  the  range  from  12  to  36  hr  was  made.  The  final  result  ob¬ 
tained  was 


A.  = 
i 

A  .  = 

i 

where  i  is  the  channel 
cutoff  in  degrees,  and 
from  the  fourth  order 


Cij(T+24)rl 


5 


0  <  T  <  12 


(3) 


£  C.jT^'1  12  STS  24 

j=l 

number,  T  is  the  local  time  in  hours,  Ais  the 
C  is  the  appropriate  coefficient  determined 
fit.  Channels  2  through  9  yielded  satisfactory 
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FIG.  17.  NINTH  ORDER  FIT  TO  THE  AVERAGE  CUTOFF  versus  LOCAL  TIME  PLOT 
REDUNDANTLY  EXTENDED  OVER  48  HOURS. 
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results  from  this  fitting  process.  The  values  of  the  coefficients  ap¬ 
pear  in  Table  II.  The  Channel  1  fit,  however,  became  physically  un¬ 
realistic  probably  because  of  the  few  data  points,  large  deviations  and 
small  diurnal  variation.  The  Channel  8  fit  is  shown  in  Figure  18. 

Several  distinct  features  show  up  in  Figure  19,  which  shows  all 
eight  fits  together.  The  local  time  plots  for  all  the  energies  seem  to 
have  maxima  around  1100  hr  local  time  and  minima  around  2300  hr. 

The  plots  seem  approximately  symmetrical  with  respect  to  the  1100- 
2300  hr  meridian  and  show  a  diurnal  variation,  which  decreases  with 
increasing  energy.  The  higher  energy  electrons  cut  off  at  lower  la¬ 
titudes  for  all  local  times;  the  spread  of  cutoffs  is  largest  around 
1100  and  smallest  around  2300  hr.  In  general,  the  cutoffs  agreed  well 
with  those  reported  by  others. 

Since  the  diurnal  variation  appeared  to  be  a  function  of  energy, 

plots  ofAAas  a  function  of  Aw  .  Alf,  .  and  A  were  run  off  for 
r  MAX  MIN'  AVG 

matched  passes  over  the  poles.  More  specifically  the  values  were 
obtained  using  only  two  data  points  taken  over  the  same  pole  on  the  same 
revolution.  This  was  designed  to  determine  whether  the  diurnal  vari¬ 
ation  was  actually  a  function  of  energy  or  a  function  of  the  magnetic 
field  configuration  where  the  energies  involved  characterize  different 
L  shells.  Williams  and  Mead  reported  a  correlation  between  the  day- 
side  cutoffs  and  the  magnitude  of  the  diurnal  variation  for  matched 
pass  data  in  1965  (Ref.  16:3021).  All  available  matched  pass  data  from 
the  present  study  was  plotted  and  each  pair  of  points  grouped 
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COEFFICIENTS  FOR  THE  EXPRESSION  OF  THE  CUTOFF  AS  A  FUNCTION  OF  LOCAL  TIME 
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into  a  2-hr  interval  according  to  the  earlier  local  time.  For  example, 
two  points  with  local  times  of  0200  and  1400  hr  would  fall  into  die 
0000-0200  interval.  The  two  local  times  sampled,  however,  were  u- 
sually  somewhat  less  than  12  hr  apart.  In  most  cases  the  two  times 
were  6  to  10.5  hr  apart. 

Many  of  the  plots  had  tco  few  points  to  indicate  any  definite 

trends,  but  examples  of  some  of  the  more  successful  plots  appear  in 

Figures  20  through  22.  The  plots  of  AA  versus  A  show  little  but 

MIN 

those  of  A  A  versus  A  show  most  points  falling  below  a  curve  with 

MAX 

a  positive  slope  and  an  X  intercept,  which  varies  inversely  as  the 

energy.  One  would  expect  points  to  fall  below  such  a  curve,  since 

the  points  represent  data  with  varying  local  time  differences,  the 

higher  points  having  differences  closer  to  the  maximum  of  12.  The 

AA  versus  Aa,.„  plots  were  similar  to  those  of  A  A  versus  A  ,  but 

AVG  MAX 

they  had  a  greater  slope  and  lower  X  intercept. 

The  results  of  these  observations  are  inconclusive.  One  would 

expect  the  positive  A  •  slope  because  AAis  dependent  directly  on 

MAX 

A  and  A  .  The  intercept  variation  may  be  due  to  the  difference 
MAX  MIN 

between  the  two  local  times  sampled  on  each  matched  pass,  which 
varies  according  to  energy  or  to  where  the  cutoff  occurred.  Since 
the  higher  energy  cutoffs  occur  at  lower  latitudes,  the  two  times  sam¬ 
pled  are  usually  farther  apart.  Thus,  one  would  expect  the  higher 


energy  curves  to  fall  above  those  of  the  lower  energies,  as  is  the  case. 
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PLOT  OF  DIURNAL  VARIATION  OF  CUTOFF  AS  A  FUNCTION  OF  THE  HIGHER 
CUTOFF  FOR  MATCHED  PASS  DATA  FROM  ALL  CHANNELS  WITH  ONE  DATA  POINT 
FALLING  IN  THE  LOCAL  TIME  INTERVAL  OF  0000-0200  HOURS 


IT  OF  DIURNAL  VARIATION  OF  CUTOFF  AS  A  FUNCTION  OF  THE  LOWER 
OFF  FOR  MATCHED  PASS  DATA  FROM  ALL  CHANNELS  WITH  ONE  DATA  POINT 
LING  IN  THE  LOCAL  TIME  INTERVAL  OF  0000-0200  HOURS 
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The  plots  would  be  expected  to  fall  somewhere  between  the 

A  and  A  results,  and,  as  far  as  the  data  indicates,  this  does 
MIN  MAX 

occur. 

The  effect  of  the  variation  of  local  time  differences  with  energy 
could  be  eliminated  by  application  of  a  selection  criterion  to  the  points 
before  plotting.  This  was  tried  in  two  ways.  Requiring  a  time  differ¬ 
ence  of  between  10  and  12  hr  reduced  the  number  of  usable  data 
points  so  drastically  that  the  plots  were  virtually  useless,  while  re¬ 
quiring  a  difference  between  9  and  12  hr  yielded  essentially  the  same 
results  as  those  obtained  with  no  selection. 

With  the  previously  obtained  expressions  for  the  cutoff  as  a 
function  of  local  time,  the  data  for  Channels  2  through  8  were  norma¬ 
lized  to  68  deg  and  replotted  as  a  function  of  the  3-hr  average  Kp  in¬ 
dices,  as  shown  in  Figure  23.  The  plots  were  made  with  varying  time 
lags.  For  example,  A  was  plotted  as  a  function  of  the  Kp  value  within 

which  the  observation  was  taken,  the  K  value  preceding  the  one  du- 

P 

ring  the  observation,  and  so  on.  The  only  trend  detected  was  that 

the  spread  of  cutoffs  was  greatest  during  times  of  higher  K  values. 

P 

To  be  sure  nothing  was  being  masked  by  the  scattering  of  the  data, 
average  values  of  the  cutoffs  were  calculated  and  plotted  over  inter¬ 
vals  of  K  .  As  shown  in  Figure  24,  nothing  new  became  apparent. 

P 

In  1968,  Williams,  Arens,  and  Lanzerotti  found  an  excellent 
correlation  between  the  appearance  of  energetic  electrons  during 
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FIG.  23.  PLOT  OF  LOCAL  TIME  CORRECTED  CUTOFF  versus  THE  3  HOUR  AVERAGE  Kp 
INDEX  FOR  CHANNEL  7  WITH  NO  TIME  LAG 


44 


-  •  r 


JF*t* 


GNE/PH/69-4 

storms  and  Dst  (Ref.  17:5673).  After  finding  no  correlation  between 

K  and  the  cutoffs.  Dst  was  tried.  Dst  is  a  measure  of  the  magnetic 
P 

field  produced  by  currents  flowing  at  great  distances  from  the  earth 
and  is  often  associated  with  ring  current  measurements.  Magnetic 
storms  produce  large  negative  Dst  values  (Ref.  10:40).  Hourly  values 
of  Dst  were  obtained  through  the  courtesy  of  Drs.  M  Suguira  and  S. 
Hendricks  enabling  plots  of  A  versus  Dst  to  be  run  off.  The  cutoffs 
were  corrected  for  local  time,  and  the  plots  were  made  for  Dst  values 
from  5  hr  before  the  cutoff  vas  observed  to  5  hr  after.  One  of  the 
plots  appears  in  Figure  25.  The  bunching  of  points  made  these  initial 
plots  difficult  to  interpret,  so  the  cutoffs  were  averaged  over  intervals 
of  Dst  and  replotted,  as  shown  in  Figure  26. 

For  Dst  values  from  -35  to  +25  gamma,  where  the  density  of 
points  is  high  enough  to  make  the  averages  meaningful,  no  trends  are 
visible,  as  can  be  seen  in  Figure  26  which  can  be  considered  repre¬ 
sentative  of  the  other  Channel  9  plots.  In  the  region  below  -35  gamma, 
where  there  are  few  points,  the  average  values  seem  to  tail  off. 

This  is  also  true  of  the  other  plots.  Figure  27  shows  the  averages  for 
all  10  of  the  correlations  for  Dst  values  from  5  hr  before  to  5  hr 
after  the  observation.  Figure  27  is  for  the  Channel  9  data  but  the 
other  channels  behaved  in  essentially  the  same  way. 
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FIG.  25.  PLOT  OF  LOCAL  Tl  ME  CORRECTED  CUTOFF  versus  DST  FOR  CHANNEL  9  WITH  NO  TIME  LAG 
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IV.  Conclusions 


Since  the  cutoffs  appesr  to  be  independent  of  altitude*  the  data 
is  apparently  not  seriously  biased  by  background  or  flux  levels.  Since 
the  flux  levels  are  functions  of  altitude  or  of  the  distance  down  the 
field  line  at  which  the  point  of  observation  is  located*  the  cutoffs 
should  show  a  dependence  on  altitude  if  the  flux  levels  are  influencing 
them  in  some  systematic  way. 

The  symmetry  with  respect  to  die  1100-2300  hr  meridian  in  die 
plots  of  A  versus  local  time  is  probably  due  to  die  angle  of  solar  wind 
interaction  with  die  magnetosphere.  In  1964*  Walters  made  a  quan¬ 
titative  estimate  of  die  angle  between  the  earth-sun  line  and  the  axis 
of  solar  wind  interaction  with  the  magnetosphere.  The  oblique  angle 
made  by  the  spiral  interplanetary  magnetic  field  vith  radially  expan¬ 
ding  solar  wind  was  estimated  to  result  in  an  easterly  deflection  of 
the  solar  wind  of  from  5  deg  to  20  deg.  The  earth’s  orbital  velocity 
relative  to  the  radial  solar  wind  direction  was  estimated  to  contribute 
an  additional  5  deg.  These  two  effects  would  cause  the  axis  of  sym¬ 
metry  of  the  magnetosphere  to  be  tilted  away  from  the  earth-sun  line 
as  if  the  solar  wind  comes  from  a  direction  10  deg  to  25  deg  west  of 
the  sun  as  shown  in  Figure  28  (Ref.  15:1769).  Thus  one  could  expect 
the  local  time  symmetry  axis  to  be  shifted  by  approximately  an  hour, 
as  is  the  case.  Strong,  however*  reports  from  Vela  satellite  obser¬ 
vations  that  at  17  Re  only  a  2  to  4  deg  east-west  tilt  of  the  magneto- 
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sphere  exists.  He  predicts  this  to  be  the  result  of  the  earth's  orbital 
velocity  abo>t  the  sun.  The  additional  tilt  from  the  interplanetary  mag¬ 
netic  field  does  not  seem  to  be  present  in  his  data  (Ref.  4:377).  Thus 
there  is  some  question  as  to  whether  part  of  the  observed  tihift  of  the 
local  time  symmetry  axis  is  due  to  the  interplanetary  field  or  to  some 
other  effect. 

No  appreciable  dawn-dusk  asymmetry  was  detectable  in  the  plots 
of  A  versus  local  time.  In  1968.  when  McDiarmid  and  Burrows  repor¬ 
ted  such  a  dawn-dusk  asymmetry  in  their  data,  they  were  dealing  with 
’Electrons  E  >  40  keV,  magnetic  local  time*  and  periods  of  little 
magnetic  activity  (Ref.  7:49).  It  is  most  probable  that  their  asym¬ 
metry  is  caused  by  electrons  with  energies  between  40  and  225  keV, 
electrons  too  low  for  this  instrument.  There  is  also  the  more  remote 
possibility  that  the  increased  scattering  of  the  points  from  variations 
in  magnetic  activity  or  the  use  of  geographic  rather  than  magnetic 
local  time  obscured  the  effect. 

The  fact  that  no  relationship  between  and  A  was  found  leads 

one  to  conclude  that  either  K  is  not  a  good  measure  of  the  distortion 

P 

v  the  geomagnetic  field  lines  or  the  3 -hr  average  value  for  Kp  is  too 

broad  in  time  to  give  a  good  correlation.  That  A  seems  to  show  the 

greatest  variation  during  times  of  high  K  incidates  that  perhaps  the 

P 

latter  is  true.  However,  the  relationship  between  A  and  K  is  still  not 

P 

clear. 
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The  tailing  off  of  the  A  versus  Dst  plots  for  large  negative  va¬ 
lues  of  Dst  suggests  that  Dst  may  be  a  good  general  index  for  the  ma* 
gnetospheric  configuration.  Hov  jver,  it  does  not  app*»«r  to  be  accurate 
as  an  exact  index,  since  correlations  with  Dst  from  several  hours  be¬ 
fore  to  several  hours  after  the  observation  yielded  essentially  the 
same  results.  Dst  is  believed  to  be  mainly  due  to  the  establishment 
of  a  ring  current  but  also  to  include  the  magnetic  field  due  to  cur¬ 
rents  on  the  Interface  between  the  magnetosphere  and  solar  plasma 
and  to  some  residue  of  the  polar  disturbances  (Ref.  10:35).  Most  suc¬ 
cessful  correlations  with  Dst  have  been  during  periods  of  high  mag¬ 
netic  activity.  In  this  study  all  levels  of  magnetic  activity  are  con¬ 
sidered.  and,  in  fact,  the  best  correlation  was  obtained  for  large  ne¬ 
gative  Dst  values.  If  the  data  were  available,  one  could  probably  ob¬ 
tain  a  better  and  more  meaningful  correlation  by  looking  at  continual 
cutoff  data  through  a  storm  as  a  function  of  Dst.  This,  however, 
would  require  nearly  complete  data  coverage. 

The  A  versus  Dst  plots  appeared  about  the  same  for  all  channels. 
If  Dst  can  be  considered  as  a  general  index  of  magnetospheric  con¬ 
figuration,  this  result  suggests  that  the  different  energy  cutoffs  are 
actually  labeling  different  L  shells. 

The  energies  examined  in  this  study  are  too  high  to  be  appre¬ 
ciably  affected  by  electric  fields  of  the  magnitude  predicted  to  be 
present  in  the  magnetosphere.  For  example,  the  Taylor  and  Hones 
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models  predicts  fields  that  would  appreciably  affect  only  particles 
with  energies  under  approximately  100  keV  (Ref.  11:3605).  Taylor 
(private  communication)  suggested  that  the  observed  energy  depen¬ 
dence  of  the  cutoffs  is  entirely  due  to  the  limiting  effect  of  rigidity  on 
the  night  side  of  the  magnetosphere.  A  particle  will  remain  trapped 
as  long  as  the  Lorentz  force  or  stabilizing  force  from  the  particle’s 
cyclotron  motion  is  greater  than  any  perturbing  force  on  the  particle. 
However,  if  the  centrifugal  force  experienced  by  the  particle  as  it 
follows  the  curved  lines  of  the  magnetosphere  exceeds  the  Lorentz 
force,  the  particle  will  be  scattered  in  pitch  angle  and  possibly  lost. 

A  condition  for  this  loss  mechanism  to  occur  is 

mv2  cos^Ot  eBv  sing  (4) 

r  c 

c 

where  m  is  the  particle’s  relativistic  mass,  v  is  its  velocity,  fl*is  its 
pitch  angle,  e  is  its  charge,  B  is  the  magnitude  of  the  field,  c  is  the 
speed  of  light,  and  r£  is  the  radius  of  curvature  of  the  field  lines.  As 
eq.  (4)  indicates,  this  mechanism  becomes  important  where  the  ra¬ 
dius  of  curvature  of  the  field  lines  becomes  small,  and  it  affects  par¬ 
ticles  with  small  equatorial  pitch  angles,  which  mirror  at  high  lati¬ 
tudes. 

In  the  magnetospheric  models  incorporating  a  current  sheet, 
the  field  lines  on  the  night  side  show  a  kink,  which  becomes  more 
severe  as  the  first  open  field  line  is  approach-  i.  This  effect  is  shown 
in  Figure  29.  The  result  is  that  the  limiting  condition  for  trapping 
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occurs  at  midnight*  where  the  curvature  of  the  field  lines  is  sharpest. 
In  applying  the  above  trapping  criterion  to  his  model*  Taylor  found 
the  theoretical  curve  shown  in  Figure  30  for  the  cutoff  as  a  function  of 
energy.  Also  shown  in  Figure  30  is  the  experimental  data  from  the 
present  study. 

The  close  correlation  between  the  slopes  of  die  observed  data 
and  predicted  curve  supports  the  conclusion  that  rigidity  is  the  mecha¬ 
nism  causing  the  observed  energy  dependence.  The  difference  in  mag¬ 
nitude  between  them  probably  exists  because  the  model  does  not  ac¬ 
curately  define  the  location  of  the  last  closed  field  line.  More  recent 
data  from  the  spectrometer  aboard  the  OVl-14  satellite  indicate  that 
the  cutoffs  for  energies  below  »80  keV  begin  occurring  at  constant  or 
lower  latitudes.  An  early  malfunction  of  the  OVl-14  has  limited  the 
number  of  orbits  of  data  available,  so  that  no  extensive  study  of  it  can 
be  made.  This  early  result,  however,  implies  that  the  E  >  40  keV 
electron  cutoffs  reported  by  previous  experimenters  are  associated 
directly  with  some  energy  around  80  keV.  The  models  derived  from 
the  E  >  40  keV  data  would  be  correspondingly  affected,  in  particular 
the  theoretical  curve  in  Figure  30.  In  any  case,  the  neutral  sheet 
models  of  the  magnetosphere  appear  valid  in  that  field  lines  are  dis¬ 
torted  in  the  predicted  manner. 

If  the  energy  dependence  is  due  to  rigidity,  the  lack  of  an  ap¬ 
preciable  dawn-dusk  asymmetry  in  the  data  suggests  that  there  is  no 
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FIG.  30.  PREDICTED  AND  OBSERVED  HIGH  LATITUDE  ELECTRON  CUTOFFS  AS  A  FUNCTION  OF 
ENERGY 
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appreciable  source  of  high  energy  electrons  in  the  sunlit  hemisphere. 
Also,  it  appears  that  energetic  electrons  do  not  drift  radially  an  ap¬ 
preciable  amount  during  one  drift  period. 

The  apparent  dependence  of  the  diurnal  variation  upon  energy 
is,  at  least  in  part,  due  to  geomagnetic  field  distortion.  Figure  31  is 
a  plot  of  the  diurnal  variation  from  the  fits  to  the  A  versus  local  time 
plots  as  a  function  of  the  day  and  night  cutoffs  for  the  lowest  eight 
energies.  Also  plotted  are  the  diurnal  variations  as  a  function  of  the 
daytime  cutoff  predicted  by  Williams  and  Mead,  using  a  40  y  finite 
current  sheet  model,  and  by  the  model  of  Taylor  and  Hones.  If  par¬ 
ticle  drift  is  independent  of  energy,  the  poor  correlation  between  the 
observed  data  points  and  the  theoretical  curves  indicates  that  these 
models  are  not  valid  in  this  critical  part  of  the  magnetosphere. 

F  igure  31  shows  that,  in  order  to  fit  the  observed  data,  the 
true  value  of  the  neutral  sheet  in  the  models  should  be  somewhere 
between  the  two  values.  However,  since  the  sizes  of  the  neutral 
sheets  are  rather  large  at  present  as  compared  to  those  actually  ob¬ 
served,  a  more  correct  model  could  probably  be  obtained  if  ring 
current  effects  were  included.  Outside  the  ring  current,  field  lines 
would  be  extended.  They  would  be  extended  to  a  greater  extent  on 
the  nightside  hemisphere,  producing  large  diurnal  variations  without 
requiring  the  high  intensity  neutral  sheets  presently  being  used. 

(Ref.  16:3027). 
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V.  Suggestions  for  Future  Investigation 

Since  this  study  did  not  reveal  any  well  defined  relationship 
between  the  high  latitude  cutoffs  and  magnetic  activity,  it  appears  that 
some  other  approaches  might  be  more  fruitful.  Perhaps  with  much 
greater  coverage,  continual  cutoffs  could  be  looked  at  as  a  function  of 
energy  through  a  magnetic  storm,  to  determine  more  exactly  what  the 
boundary  is  d*ung  during  disturbed  times.  Also  one  could  investigate 
the  effect  of  the  ring  current  more  thoroughly  by  looking  at  this  con¬ 
tinual  data  along  with  Dst  and  perhaps  even  the  rate  of  change  of  Dst 
with  time. 

The  mechanism  thought  by  Taylor  to  be  responsible  for  the  ob¬ 
served  energy  dependence  of  the  cutoff  would  result  in  cutoffs  that 
are  independent  of  the  maximum  outer  zone  flux.  Thus  one  could  test 
the  theory  by  determining  whether  a  correlation  does  exist  between 
the  cutoff  and  the  maximum  flux  value.  The  theory  could  also  be 
tested  by  determining  the  altitude  dependence  of  the  cutoff  along  a 
field  line  in  the  area  where  the  cutoff  is  observed  to  occur.  Since  the 
loss  mechanism  is  dependent  on  pitch  angle,  the  flux  in  the  equatorial 
plane  should  be  well  above  background  on  a  field  line  where  a  high 
latitude  cutoff  occurs. 

The  energy  dependence  of  the  diurnal  variations  in  cutoffs  left 
some  question  as  to  whether  or  not  particles  follow  drift  paths  inde¬ 
pendent  of  energy.  Construction  and  examination  of  isointensity  plots 
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should  aid  in  answering  this  question. 

Behavior  of  the  cutoffs  as  a  function  of  energy  for  electrons  of 
energy  E  <  225  keV  has  not  been  examined  in  detail.  The  OVl-19 
should  provide  the  necessary  data  for  such  a  study.  Low  energy  data 
is  essential  for  determining  electric  field  configuration  in  the  magne¬ 
tosphere  as  well  as  particle  source,  loss,  diffusion,  and  acceleration 
mechanisms. 

The  high  latitude  spikes  noted  in  reading  off  the  cutoffs  were 
not  investigated  in  this  study.  They  are  real,  however,  and  should  be 
investigated  in  the  future. 

The  electron  flux  data  used  in  this  study  now  exists  in  the  form 
of  plots  of  flux  versus  L.  This  introduces  the  possibility  of  studying 
the  cutoff  mechanism  more  closely.  Such  parameters  as  the  slope  of 
the  flux  versus  L  curve  might  be  worthy  of  investigation. 
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Appendix  A:  The  Instrument 

The  instrument  used  was  the  magnetic  spectrometer  shown  in 
Figure  8.  Electrons  enter  the  indicated  aperture  and  are  bent  by  the 
magnetic  field  according  to  their  energy.  An  electron  passing  perpen¬ 
dicularly  through  a  B  field  will  circle  about  a  field  line  according  to 
the  relation* 

»  .22.  -S-/z  m£  +E?  V/2  (5) 

c  eB  eB\  c  )  '  ’ 

where  is  ihe  particle’s  radius  curvature,  m  is  its  relativistic 
mass,  v  is  its  velocity,  perpendicular  to  ihe  field,  e  is  its  electronic 
charge,  E  its  kinetic  energy,  c  is  the  speed  of  light  and  B  is  the  value 
of  the  field.  Thus  particles  of  different  energies  and  masses  will  have 
different  radii  of  curvature. 

The  magnetic  spectrometer  used  to  acquire  the  data  was  of  the 

180-deg  type  with  nine  solid  state  detectors  located  at  the  primary 

focus.  The  detectors  were  1.0cm  x  1.5cm  x  1000//  thick  lithium 

drifted  diodes  with  a  0.95cm  x  1.37cm  unshielded  sensitive  area.  Each 

detector  had  its  own  amplifier,  twc  level  pulse  height  analyser  and  two 

log  pulserate-to-analog  converters  with  useful  ranges  from  1  to  10 

5 

and  from  300  to  3.0  x  10  pulses  per  second.  All  lower  thresholds 
were  set  at  1 50  keV  and  all  but  one  upper  threshold  were  set  at  2.65 
MeV  in  order  to  get  uniform  bremsstrahlung  and  penetrating  proton 
background  in  each  channel.  The  upper  threshold  in  Channel  1  was 
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set  at  2.9  MeV.  The  energy  deposits  above  2.65  MeV  in  Channel  9 
were  monitored  once  each  second  to  determine  energetic  proton  con* 
tributions.  Minimum  energy  for  direct  proton  penetration  is  105  MeV. 

Internal  baffling  around  the  edges  of  the  uniform  field  volume 
on  the  pole  plates  acts  as  a  collimator  to  reduce  scattering  within 
the  instrument.  During  calibration  no  scattered  electrons  were  ob¬ 
served,  though  bremsstrahlung  sensitivity  is  sufficiently  low  for  the 
spectrums  being  investigated  and  no  corrections  to  the  data  were 
necessary. 

A  4-MeV  Van  de  Graaff  accelerator  was  used  to  check  energy 
levels,  angular  acceptance,  and  bremsstrahlung  sensitivity.  Detector 
efficiencies  were  determined  using  electrons  from  Ru*°^  source  in 
a  laboratory  180-deg  magnetic  spectrometer.  Protons  up  to  150  MeV 
were  used  to  determine  the  instrument  response  to  penetrating  ener¬ 
getic  protons.  Final  calibration  of  energy  levels  and  geometric  fac¬ 
tors  was  done  by  trajectory  analysis  on  a  computer  using  the  mea¬ 
sured  field  and  fringing  field  geometries  and  the  physical  dimensions 
of  the  instrument.  Actual  particles  were  used  to  verify  the  compu¬ 
ter-generated  energy-angle  response  curves.  The  various  charac¬ 
teristics  of  the  nine  channels  appear  in  Table  m.  All  internal  vol¬ 
tages,  detector  currents,  and  internal  temperatures  are  sampled  once 
each  eight  seconds  primarily  to  determine  the  state  of  the  detectors. 
No  changes  have  yet  been  observed,  probably  due  to  the  high  thres- 
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TABLE  III.  Characteristics  of  Detectors 


Detector 

CENTER  E(MeV) 

AE  (MeV) 

Geometric  Factor 

(cm2  ster) 

1 

2.31 

0.325 

2.27  x  10“2 

2 

2.04 

0.320 

2.48  x  10~2 

3 

1.77 

0.310 

2.77  x  10~2 

4 

1.49 

0.322 

3.10  x  10“2 

5 

1.225 

0.300 

4.17  x  10-2 

6 

0.957 

0.285 

4.67  x  10'2 

7 

0.712 

0.275 

5.46  x  10’2 

8 

0.475 

0.250 

7.30  x  10~2 

9 

/  0.300 

\  >105  protons 

0.150 

4.63  x  10"2 

(Ref.  14:3-7) 
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hold  for  proton  penetration.  The  instrument  is  still  operating  as  it 
did  in  the  laboratory  during  calibration  except  for  Channel  2  which  be** 
came  noisy  after  18  months  in  orbit.  The  data  used  in  the  study  was 
taken  well  before  this  malfunction  occurred. 

The  spacecraft  was  spin- stabilized  at  8.8  rpm.  The  particle  de¬ 
tector  was  mounted  with  its  field  of  view  perpendicular  to  the  spin 
axis  and  so  spin  provided  a  scanning  function.  An  on-board  triaxial 
fluxgate  magnetometer  provided  local  magnetic  field  data  which  was 
used  to  determine  pitch  angles  for  the  directional  data.  A  1 50-minute- 
capacity  tape  recorder  permitted  data  to  be  acquired  over  an  entire 
orbit  (Ref.  14:1254-1256). 
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Appendix  B:  Error  Propagation 

The  errors  associated  with  determining  the  universal  time  at 
which  the  cutoff  occurs  are  primarily  of  the  random  type.  The  great¬ 
est  source  is  the  error  of  judgment  associated  with  finding  the  cut¬ 
offs.  Since  the  amount  of  error  arising  from  this  source  is  largely 
dependent  upon  the  sharpness  of  the  cutoff  and  the  density  of  data 
points  along  the  curve,  it  varies  according  to  each  reading.  To  ob¬ 
tain  an  average  standard  deviation  for  all  the  data  would  involve  the 
difficult  time-consuming  task  of  reading  off  all  the  cutoffs  several 
times.  Some  of  the  more  questionable  cutoffs,  however,  were  picked 
off  and  the  consistency  of  the  readings  checked.  This  showed  that  the 
cutoffs  as  deiined  had  at  most  an  error  of  *30  seconds  due  to  judg¬ 
ment.  The  corresponding  error  in  terms  of  invariant  latitude  depends 
oil  how  fast  tiie  satellite  is  moving  and  where  it  is  at  the  time  the  cut¬ 
off  is  determined.  On  an  apogee  pass  over  a  pole  the  30  seconds 
would  correspond  to  0.5  deg  invariant  latitude  while  on  a  perigee  pass 
it  could  be  as  high  as  2.2  deg. 

Another  source  of  random  error  arises  from  the  ephemeris. 

The  ephemeris  for  some  of  tiu<  data  may  be  poor  and  introduces  a 
possible  error  of  up  to  10  seconds  in  universal  time  and  a  corres¬ 
ponding  error  in  the  cutoff  latitudes. 

The  only  known  possible  source  of  systematic  error  is  biasing 
of  the  data  by  the  observer  because  of  the  cutoff  shapes.  It  is 
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possible  that  the  shape  of  the  curve  around  the  cutoff  as  influenced 
by  the  maximum  flux  and  slope  influenced  the  reading  in  some  way. 
We  know  of  no  such  error  but  do  introduce  the  possibility.  In  any 
event  its  contribution  is  small  in  comparison  to  the  random  sources. 
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